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ABSTRACT 

We use the large catalogues of haloes available for the Millennium Simulation to test whether 
recently merged haloes exhibit stronger large-scale clustering than other haloes of the same 
mass. This effect could help to understand the very strong clustering of quasars at high red- 
shift. However, we find no statistically significant excess bias for recently merged haloes over 
the redshift range 2 ^ z ^ 5, with the most massive haloes showing an excess of at most ^5%. 
We also consider galaxies extracted from a semianalytic model built on the Millennium Simu- 
lation. At fixed stellar mass, we find an excess bias of 20 — 30% for recently merged objects, 
decreasing with increasing stellar mass. The fact that recently-merged galaxies are found in 
systematically more massive haloes than other galaxies of the same stellar mass accounts for 
about half of this signal, and perhaps more for high-mass galaxies. The weak merger bias of 
massive systems suggests that objects of merger-driven nature, such as quasars, do not cluster 
significantly differently than other objects of the same characteristic mass. We discuss the im- 
plications of these results for the interpretation of clustering data with respect to quasar duty 
cycles, visibility times, and evolution in the black hole-host mass relation. 

Key words: cosmology: theory - cosmology: dark matter - galaxies: formation - galaxies: 
high-redshift - quasars: general 



1 INTRODUCTION 

In the last decade, much theoretical work has tried to constrain 
the cosmological evolution of supermassive black holes (BHs) 
by simultaneously interpreting the statistics of quasars/BHs and 
their clustering as a function of redshift and luminosity (e.g.. 



thousands of quasars, has allowed a detailed investigation of the 
clustering properties of accreting BH s from the local universe up 
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iKauffmann & Haehnell' '2002 d: IWvithe & Loebl |200^: iLidz et al.l 
2006; Hopkins et al. 2007; Th acker et al.ll2009l ; lBonoli et al.ll2009l ; 
Ishankar et al.ll2008L ,2009). In fact, if quasars are hosted by haloes 
whose bias is only mass-dependent, clustering measurements can 
be used to infer the mass Mnaio of the host dark matter halo, which 
in turn provides host number densities, quasar duty cycles (here 
defined as the ratio between quasar and halo number densities) and 
scatter in t he relation between quasar luminosity L and h alo mass 



I Cole & Kaiser. .19891; iHaehnelt et alj [1998 ; Martini & Weinberg! 



I2OOII; iHaiman & Huill200l[r Measuring a high bias implies high 
halo masses, low host number densities, high duty cycles, and vice 
versa. At fixed duty cycle, increasing the scatter in the mean L- 
Mnaio relation implies increasing the contribution of less massive 
and less biased haloes to the same luminosity bin, thus lowering the 
overall bias. 

The advent of wide-field surveys like the Sloan Digital 
Sky Surve y (SD SS) and the 2dF q uasi-stellar object (2dFQS0) 
( I York et alJi200Ql ; iCroom et al.ll2004l) survey, with the detection of 



E-mail:bonoli@mpa-garching. mpg.de 



hosting quasars are typical in the way they trace the dark mat- 
ter density field, these studies concluded that quasars reside at 
all times in a relatively narrow range of halo masses, MHaio~ 
3x 1012-1013/j-1Mr. 



Interestingly, the very high clustering amplitude of lumi nous 
quasars at z > 3 measured with the SDSS dShen et alj|2007l) has 
posed some nagging theoretical problems for the simultaneous in- 
terpretation of the clustering and the luminosity function at these 
epochs. The high clustering appears to imply that the quasars live 
in very massive haloes. But the extreme rareness of such haloes 
is difficult to reconcile with the observed quasar luminosity func- 
tion, especially at z ~ 4, unless a high quasar duty cycle and a 
very low scatter in the L-Mjjaio relation are assumed jWhite et al.l 
2008). Moreover, matching the high z > 3 — 4 quasar emissivity 
to the low number density of hosts constrains the ratio of the ra- 
diative efficiency of accretion e to the Eddington ratio X to be 
e>() .7A,/(l + 0.7X), implying e > 0. 17 for X > 0.25 jShankar et all 
|2008|) . which are rather extreme values. However, these conclu- 
sions can be relaxed if qua sar hosts cluster more s trongly than typi- 
cal haloes of similar mass jWvithe & Lo"ebll2009l) . This would then 
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imply that quasars live in less massive but more numerous haloes, 
allowing for lower duty cycles and less extreme values for £. 

Several analytical and numerical studies have investi- 
gated whether haloes of similar mass have different large- 
scale clustering properties, depending, in a non-trivial way. 
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20081). In particular. 



Wvithe&Loebl ( l2009l) suggest that the possible merger-driven na 



ture of quasars might cause an excess bias, if the large-scale clus- 
tering of recently-merged haloes is higher than expected for typi- 
cal haloes of the same mass ("merger bias"). This suggestio n was 
based on the model by iFurlanetto & Kamionkowskil j2006l) . who 
analytically calculated that close merging pairs might possess a 
merger bias of a factor of > 1.5. 

Recent work has indeed shown that clustering strength de- 
pends not only on halo mass, but also on other parameters. 
iGao et al. I ( l2005h found that later forming haloes with mass M < 
M* are less cl ustered than tjfpical h aloes of the same mass ("as- 
sembly bias"). IWechsler et alj |l2006) extended this study to show 
that less concentrated haloes more massive than the non-li near mass 
scale are instead more biased than average. lLi et al. liooi) explored 
various definitions of halo formation time, and concluded that the 
dependence of clustering on halo history depends strongly on the 
precise aspect of the history that is probed: while they confirmed 
previous results on assembly bias, they did not find any dependence 
of clustering on the time of the last major merger. Other numerical 
work that specifically looked at merger bias found inconclusive re- 
sults, probably due to the different ranges of masses, redshi fts and 
scales used and the poor statistics available jOottlober et al ]|2002l; 
IPercival et all |2003| ; IScannapieco & Thacked l2003h . IWetzel et al.l 
( l2007h used a large dark matter simulation to study the clustering 
of very massive haloes (Mnaio > 5.0 x lO'-^/;^' Mq ), and found 
that, at redshift z < 1, merger remnants show an excess bias of 
~ 5 - 10%. 

In this work we make use of the l arge, publicly availab le cat- 
alogues of the Millennium Simulation jSpringel et alj|200^ to ex- 
plore for a wider range of masses the level of excess bias for high 
redshift merger remnants. In Section[2]we briefly describe the simu- 
lation and how we identify recent mergers both of haloes and galax- 
ies. In Section[3]we show results for the merger bias both of haloes 
and of galaxies, and in Section|4]we discuss the implications of our 
results for the clustering of quasars. A summary and our conclu- 
sions are presented in Section|5] 



2 IDENTIFYING MERGING HALOES AND GALAXIES 
IN THE MILLENNIUM SIMULATION 

2.1 The Millennium Simulation and its galaxy population 

The Millennium Simulation jSpringel et alj |200^ is an N-body 
simulation which follows the cosmological evolution of 2160"' ~ 
10^" dark matter particles, each with mass ~ 8.6 x lO^/i^'Me, in 
a periodic box of 500/i^'Mpc on a side. The cosmological parame- 
ters used in the simulation are consistent with the WMAPl &. 2dF- 
GRS 'concordance' ACDM framework: fl„, = 0.25, = 0.75, 
08 = 0.9, Hubble parameter h = //o/100 kms~'Mp c~' = 0.73 and 
primordial spectral index n=l jSpergel et al.l2003l) . In the present 
work we focus on the clustering of galaxies and haloes from z = 2 to 
z = 5, where the time between two simulation outputs varies from 



approximately 200 Myr to 100 Myr. This time resolution is good 
enough to capture merger events reliably and in these time inter- 
vals any change in the large-scale distribution of merger remnants 
is negligible. 

Detailed merger trees were constructed for the simulation by 
identifying haloes and subhaloes with, respectively, a friends-of- 
friends (FOF) grou p-finder a nd an extended version of the SUB- 
FIND algorithm ( ISpringel et al.ll200lh : particles are included in the 
same FOF group if their mutual separation is less than 0.2 of the 
mean particle separation. The SUBFIND algorithm then identifies 
locally overdense and self-bound particle structures within FOF 
groups to isolate bound subhaloes (which are required to contain 
a minimum of 20 particles). For further details on the Millennium 
Simulation and the t ree building procedure we refer the reader to 
ISpringel et al] ( l2005h . 

The formation and evolution of galaxies has been followed in 
a post-processing simulation which uses the dark matter merger 
trees as basic input combined with analytical treatments of the 
most important baryonic p hysics in g alaxy for mation, includin! 
the growth of central BH s jCroton et al. 2006; Bower et al1l200i 
iDe Lucia & Blaizo^l2007h . This has produced remarkably success 
ful galaxy formation models which reproduce a large set of obser- 
vational findings about the local and high redshift galaxy popula- 
tions with good accuracy. While not perfect, this match justifies 
substantial trust in the basic paradigm of hierarchical galaxy for- 
mation in CDM cosmologies, and motivates detailed studies of the 
merger and clustering statistics using the Millennium Simulation. 
Below we describe our definition of major mergers both for dark 
matter haloes and galaxies, which lies at the heart of our investiga- 
tion of the merger bias phenomenon. 

2.2 Halo mergers 

We note that different definitions of halo formation time have led 
to somewhat different quantitative conclusions regarding the effect 
of assembly history on the large-sca l e clustering o f haloes (e.g., 
iGao et al.ll200^ ; IWechsler et al.ll2006l ; lLi et alj|2008h . Here we are 
interested in the possible bias caused by recent merger activity, 
which might induce a non-trivial relation between the clustering 
of dark matter haloes and objects whose formation is triggered by 
mergers (such as quasars). We are therefore not interested in track- 
ing the full mass accretion history of dark matter haloes, but rather 
want to focus on the violent major merger events that are thought 
to trigger efficient BH accretion and starburst activity. 

In the present work, we consider as major mergers those 
events in which two separate haloes with comparable masses en- 
counter each other for the first time, that is, when they join the 
same FOF group. At a given time Zn, we define as recently merged 
haloes those that, at the previous snapshot z„-i, have two or more 
progenitors belonging to separate FOF groups whose mass (de- 
fined through the number of particles of the FOF group) was 
> 20% that of the descendant (corresponding to a minimum ratio 
OTsat : "iccn tral = 1:4). This definition of a m ajor merger is similar to 
the one of IScannapieco & Thacked ( |2003|) . who defined as merger 
remnants haloes that, within the time interval of a single snapshot, 
accreted at least 20% of their final mass. These authors, however, 
also considered haloes that experienced considerable smooth mass 
accretion, whereas we strictly require the merger remnant to be the 
product of the en counter of two suffic iently massive FOF progeni- 
tors. We note that I Wetzel et al.l ( |2007|) used different definitions of 
a major merger, finding their results to be independent of the exact 
definition, as long as a substantial amount of mass is accreted in 
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a relatively short time. We therefore also expect our results to be 
robust with respect to the specific value of the mass ratio threshold 
adopted here. 

In our definition of major mergers we also include encoun- 
ters of groups that at some later time might split again. This can 
occasionally happen since the FOF algorithm sometimes links two 
haloes that are just passing close to each other but that in the fu- 
ture will (at least temporarily) separate again. To check whether 
this might impact our overall results, we also u sed the merger tree s 
extracted from the Millennium Simulation bv lGenel et al.l JzOOSi) . 
who carefully excluded all mergers of FOF groups containing sub- 
haloes that at a future time wil l belong to two different FOF groups. 
Moreover, iGenel et alj j2008l) define as halo mass the sum of just 
the gravitationally bound particles. We ch ecked, however, tha t our 
results do not change when switching to the lOenefet ai]( l2008h halo 
trees. The differences from our reference catalogues affect the halo 
population only at very low redshifts and at low halo masses, much 
below the ranges of interest here. 

We have also checked that our definition of halo mass based 
on the number of linked particles instead of a spherical overdensity 
mass estimate does not affect our result. As an alternative to the 
FOF group masses, we used as group masses the mass within the 
radius that encloses a mean overdensity of 200 times the critical 
density, or the mass within the radius where the overdensity is that 
expected for virialization in the generalized top-hat collapse model 
for our cosmology. However, we found that this did not lead to 
any significant differences in the large-scale clustering properties 
of haloes as a function of mass. 



2.3 Galaxy mergers 

In the Millennium Simulation, the orbits of dark matter subhaloes 
are followed until tidal truncation and stripping due to encoun- 
ters with larger objects cause them to fall below the simulation 
resolution limit (20 particles, equivalent to a mass of --^ 1.7 x 
lO'^/i^' Mq). Galaxies follow the orbits of their host subhalo un- 
til this point, and then their remaining survival time as satellite 
galaxies is estimated using their current orbit and the dynami- 
cal friction formula ofisinnev & Tremaind ( Il987h . calibrated as in 
iDe Lucia & Blaizo^ j2007h . At the end of this interval, a satellite 
galaxy is assumed to merge with the central galaxy of the host dark 
matter halo, which can either be a subh alo or, more frequen tly, the 
main halo of the associated FOF group l lAngulo et alj200§) . 

In the event of a minor galaxy merger, the cold gas of the 
satellite galaxy is transferred to the disc component of the cen- 
tral galaxy together with the stars produced in a starburst (as de- 
scribed below); moreover, the bulge of the central galaxy grows by 
incorporating all the stars of the satellite. If instead a major galaxy 
merger has occurred, the discs of both progenitors are destroyed 
and all stars in the merger remnant are gathered into the spheroidal 
bulge component. In the galaxy formation model studied here, 
the starbursts induced by galaxy mergers are de scribed using the 
"collis ional starburst" prescription introduced bv lSomerville et au 
( I2OOII) : the fraction eburat of cold gas which is converted into stars 
in the merger remnant is given by: eburst = Pburst {msm/ mcentral ) . 
where otburst = 0-7 and Pb nrst = 0.56 , chose n to provide a good fit 
to the numerical results of 



CoxetalJ ( l2004b . 



We define as major merger remnants those galaxies that have, 
in the immediately preceding simulation output, two progenitors 
with stellar masses larger than 20% of the stellar component of the 
descendant (as for the FOF haloes, this imposes a minimum mass 
ratio rasat : "Jcentral = 1:4). Note that this definition is close, but not 



identical, to the distinction between minor/major mergers adopted 
in the underlying galaxy formation model. 



3 RESULTS 

3.1 Clustering analysis and the excess bias F 

We use the standard definition of the two-point spatial correlation 
function as the excess probability for finding a pair of objects at a 
distan ce r, each in the volume elements dVi and dV2 (e.g., (Peacock 
Il999h : 



dP 



■-n^ll- 



-Ur)]dVidV2, 



(1) 



where n is the average number density of the set of objects under 
consideration. 

The bias between two classes of objects (e.g., haloes and dark 
matter) is defined as the square-root of the ratio of the correspond- 
ing two-point correlation functions: 



"H.DM 



^dm('') 



(2) 



Since the number density of merging objects at a given snap- 
shot is too low for a statistically significant auto-correlation study 
(see section|4l(, we adopt a cross-correlation analysis instead. In this 
case the bias is given by 



I 



(3) 



where i'R.DM('') is the bias (relative to the dark matter) of the pop- 
ulation we are using as reference in our cross-correlation analysis, 
and t,f{,R{r) is the cross-correlation function between the haloes 
and the reference population. By definition, the bias is a function 
of scale. However, the scale dependence becomes weak or even 
vanishes at large scales. Since we are here interested in the behav- 
ior of the merger bias at very large scales, we estimate the bias 
on these scales by finding the best constant value over the range 
5 < r < 25 /i^' Mpc. This adds robustness to our results by reduc- 
ing noise from counting statistics. 

We can define the merger bias as the excess in the clustering 
of merger remnants at large scales with respect to the global popu- 
lation of objects selected with similar properties: 



F{r) = ^M.R{r)/^H.Ri^) 



(4) 



where ^m.r is the cross-correlation between merger remnants and 
the reference sample, and ^h,r is the cross-correlation between the 
global population and the reference sample. 

We estimate errors for our measurements using the bootstrap 
method, generating for each sample 50 bootstrapped samples of 
the same size, drawn at random from the parent sample and al- 
lowing for repetitions (the error estimates converge already when 
using just a few dozen bootstrap samples). For each bootstrap sam- 
ple, we calculate the correlation functions, the bias and the excess 
bias. The standard deviatio n among these quanti ties is then taken 
as error estimate. Recently, iNorberg et alj 1 I2OO8I) pointed out that 
the variance on the two-point correlation function is overestimated 
when calculated with bootstrap techniques. Keeping this in mind, 
we deliberately choose the bootstrap method in order to be conser- 
vative in our error estimates. Another option would have been to 
estimate errors by subdividing the whole Millennium volume into 
subvolumes (for example eight octants) and then to calculate the 
variance of the ^(r) measured within individual subvolumes. This 
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Figure 1. Upper panels: Examples of the two-point cross-coiTelation function at z = 2 and z = 4 (left and right, respectively) between the reference sample 
and the haloes with mass in the range 2.0 < Mnaio < 4.0 x lO'- M,., (blue-dotted lines), and between the reference sample and the sub-sample of recently 
merged haloes (red lines). The auto-correlation of the underlying dark matter is shown as dot-dashed line. Middle panels: Bias as a function of scale for all the 
haloes in the mass bin (blue triangles), for the corresponding merger remnants (red bow-ties) and for the reference sample (black squares). The horizontal lines 
indicate the fit to the points, over the range indicated by the vertical dotted lines.. Lower panels: Excess bias F for the merger remnants relative to the whole 
galaxy population as a function of scale. The horizontal dashed line indicates F = I. We refer the reader to the text for details of the calculation of errors and 
the fitting procedure. 
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Figure 2. Bias and F parameter as a function of redshift, from the best fit 
obtained for the halo samples shown in Figure [T] 



method becomes inaccurate at large scales (few tens of Mpc) due 
to the smaller volume probed by each subvolume. 



3.2 The merger bias for DM haloes 

In our study of the merger bias for haloes we proceed as follows: 

• We take all FOF haloes with mass in the range 5 x lO'^ < 
Meaio < l-6x lO'-' /!^' Me. For the redshifts analyzed in this work, 
this mass range is well above the collapsing mass M*, defined by: 
ct(M*) = 1.69 (at z = 2, M,, ~ 1.3 X 10'° /j^'Mq). This entire sam- 
ple is used as reference sample for the cross-correlation analysis. It 
is large enough that the error on its auto-correlation can be safely 
neglected with respect to other sources of error in the b and F 
parameters (it is composed of ~ 3.5 x 10'* haloes at z = 5 up to 
~5.5x 10^ atz = 2). 

• We subdivide this sample into five mass-bins, with constant 
logarithmic spacing AlogMnajo = 0.3 (a factor of two in mass). 
We will refer to these five samples as Hi . 

• We then checked which haloes in each of the bins of had a 
recent major merger, as described in section [2!2l The subsamples of 
recently-merged objects are denoted by M, . The fraction of merger 
remnants is ~ 10% at z = 2, and increases to 15 — 20% at z = 5. 
The mass bins are narrow enough that, within each bin, the merger 
remnants and the parent population have effectively the same dis- 
tribution of masses. 

For the bootstrap error calculation, we created 50 samples from 
each of the //, halo samples, and from these new samples we ex- 
tracted the corresponding catalogues of the recently-merged haloes. 
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Figure 3. Excess bias for DM haloes in separate mass bins, as indicated in 
each panel in units of lO'^/i^' M,;). 



If ^Mj.fi is the cross correlation between the merger remnants 
and the reference sample, and ^h,,« is the cross correlation between 
all the haloes in the bin and the reference sample, the excess bias 
parameter is given by F = ^m, ,r/^h, .«- One of our principal aims 
is to quantify how much F deviates from unity. 

In Figure [T] we show at two different redshifts an exam- 
ple of the two-point cross-correlation function, the bias and the 
F parameter for haloes with masses in the range 2.0 < Mnaio < 
4.0 X lO'^/;^' Mp. In the top panels, the red and blue curves are 
the cross-correlation functions between the reference sample and 
the merger and parent halo samples, respectively. The error bars 
show the 1 a dispersion of the bootstrap samples. The correlation 
power drops at small scales as expected for FOF haloes (by def- 
inition, two haloes cannot be closer than twice their virial radius, 
hence the "1— halo" term, i.e., the contribution to the correlation 
function from subhaloes within the virial radius, is missing). In the 
middle panels, we show at each scale the bias of the merger sam- 
ple, of the parent sample and of the reference sample (red bow-ties, 
blue triangles and black squares, respectively); in the lower panels 
the excess bias F = ^m^.r/^H/.r is also shown as a function of scale. 
The errors on each point on b (and F) are from the 1 a dispersion 
of the bias (F) calculated for each bootstrap sample. Both for h and 
F the horizontal lines show the best constant fits to the points in the 
range 5 < r < 25 /j^'Mpc. 



The resulting fits for the bias and F as a function of redshift 
are shown in Figure |2] The errors on these fits are given by the 
1 a dispersion of the fits calculated for each bootstrap sample. The 
excess bias F corresponding to each halo mass bin considered is 
shown in Figure [3] If at a given snapshot there are less than 10 
mergers, we do not show the result, since the corresponding cross- 
correlation function would be too noisy. That is why for the higher 
mass-bins (lower panels) the results are not shown at all redshifts. 

In these results, we do not find any statistically significant 
merger bias, over the full redshift range probed in our analysis. 
Only in the most massive bins (lower panels), we see a possible 
signal at the < 5% level for the s mallest redshift. We stress that 
switching to the lGenel et al.l ( l2008h catalogs or changing our mass 
definition does not alter the results presented in the Figures dis- 
cussed above. 



3.3 The merger bias for galaxies 

We investigated the merger bias of galaxies with a procedure simi- 
lar to the one adopted above for dark matter haloes. All the galaxies 
with stellar mass in the range 5 x 10^ < A/star < 1.6 x lO" /i^'M© 
have been divided into five mass bins, G/, and from each bin we ex- 
tracted subsamples of recently-merged galaxies Mj, obtained as de- 
scribed in section [23l We use as reference sample the entire galaxy 
population in this range (5 x 10^ < Mstar < 1.6 x lO" /i^'Mq, 
which is composed of ~ 10^ galaxies at z = 5 up to --^ 1.4 x 10* 
at z = 2). This sample is again large enough that the error on its 
auto-correlation can be safely neglected. 

In Figure |4l we show an example of the two-point cross- 
correlation function for the intermediate mass bin. We refer to the 
description of Figure [T] for details on the derivation of the bias b 
and the excess bias F. Unlike for FOF groups, where it only makes 
sense to consider the clustering properties on large scales, for the 
galaxies we can compute the correlation function down to very 
small scales ~ 0.01 /i^'Mpc, allowing for a rather accurate descrip- 
tion of the 1-halo term as well, at least at z < 4. We find that while 
F at large scales is approximately constant, it steadily increases at 
the smallest scales probed by our study. The detection of a steady 
increase of the excess bias F with decreasing scale might be of po- 
tential interest. A comparison between model predictions and the 
observed small-scale clustering of quasars at different redshift and 
luminosity thresholds could, in fact, provide insights on the merger- 
driven nature of quasars, and we postpone a detailed analysis of this 
subject to future work. 

The excess bias F fitted on scales larger than 5 /j^'Mpc is 
plotted as a function of stellar mass and redshift in Figure [5] (solid 
lines). Excess bias of size 20 — 30% (F ~ 1.2 — 1.3) is clearly 
present for all mass bins and at all redshifts, despite the large er- 
ror bars at the highest redshifts. In essence, we find that, at fixed 
stellar mass, recently merged galaxies are more strongly clustered 
on large scales. 

We also examined the mass-distribution of the dark matter 
subhaloes hosting the galaxies considered in this analysis. For 
galaxies that sit in the main halo of a FOF group, this mass is given 
by the virial mass of the group (defined as the mass within the ra- 
dius that encloses a mean overdensity of 200 times the critical den- 
sity of the simulation), whereas for galaxies that are located in sub- 
structures, the parent halo mass is defined simply as the number of 
particles bound to the substructure (as determined by the SUBFIND 
algorithm). We found that, for each galaxy bin, the distribution of 
masses of the host subhaloes is typically log-normal, and peaks at 
systematically higher subhalo masses for recently-merged galax- 
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Figure 4. Upper panels: Examples of the two-point cross-con'elation function aX z = 2 and z = 4 (as labeled) between the reference sample and galaxies with 
stellar mass in the range 2.0 < Mstar < 4.0 x lO'" /i^'M ;, (blue-dotted lines), and between the reference sample and the sub-sample of similar mass merger 
remnants (red lines). The auto-correlation of the underlying dark matter is shown as dot-dashed line. Middle panels: Bias as a function of scale for all the 
galaxies in the mass bin (blue triangles), for the corresponding merged galaxies (red bow-ties), and for the reference sample (black squares). The horizontal 
lines show fits to the points, over the range indicated by the vertical dotted lines. Lower panels: Excess bias F for the merged galaxies relative to the whole 
galaxy population as a function of scale. The horizontal dashed line indicates an absence of excess bias (i.e., F = I ). We refer the reader to the text for details 
of the calculation of en'ors and the fitting procedure. 



ies. The median host subhalo masses for the stellar-mass bins are 
shown as a function of redshift in Figure |6] 

This raises the natural question of whether the excess bias de- 
tected for galaxies could be due simply to this offset in the typical 
mass of the host subhalo population. To test this idea, we generated 
for each galaxy bin a new parent galaxy population with the same 
distributions of stellar mass and host subhalo mass. The excess bias 
between this "corrected" galaxy population and the corresponding 
recently-merged galaxies is shown in Figure [5] as a dashed line. 
This exercise significantly decreases the excess bias signal, and no 
clear dependence on stellar mass or redshift remains. Nevertheless, 
a statistically significant excess bias (at the '--^ 10 — 20% level) still 
seems to be present, especially for the lower stellar-mass bins. 

In summary, while for FOF dark matter haloes we did not 
find any statistically significant merger bias, for galaxies a signal is 
present at a level of ~ 10 — 20% for the smallest systems. However, 
when we restrict ourselves to galaxies at the center of FOF groups 
(~ 75 - 85% of galaxies at z = 2 and ~ 95 - 98% at z = 5, de- 
pending on stellar mass), the excess bias approaches that obtained 
for dark matter haloes alone ( section |T2t . The differing results ob- 
tained for haloes and the galaxy population must therefore be due 
to the physics of the merger of galaxies, which goes beyond that 
of halo merging. This is still a topic of active research in galaxy 
formation modelling. 



4 IMPLICATIONS FOR THE CLUSTERING OF 
QUASARS 

The la rge clustering amplitude of quasars observed bv IShen et al.l 
l l2007l) at high redshift appears to suggest that these objects live in 
very massive haloes. In Figure|7] the bias associated with FOF halo 
merger events for different mass ranges and at different redshifts, is 
compared to the observed quasar biafl as calculated bv lShen et al.l 
(2009). The high observed clustering is compatible with the clus- 
tering associated with the most massive DM haloes, which, at least 
up to z ~ 4, we find to b e in better agreement with the analytica l 
predictions of ljind ( Il998b . rather than those of lSheth et al] ( 1200 1[) . 
though still somewhat higher at the highest redshifts. 

As discussed in Section [T] the high observed clustering sig- 
nal forces quasar models to adopt extreme values for some of the 
relevant parameters, such as assuming very low scatter in the L- 
MHain relation, high duty cycles, and high radiative efficiencies 
JWhite et al]|2008l : Ishankar et al.ll2008 '). However, an excess bias 
applying specifically to quasar hosts compared with random haloes 
of the same mass might reduce the need for such strong assump- 

' To correct for the diffe rent cosmo l ogies u sed, the large scale quasar 
bias measurements from IShen et al] j2009l) have been multiplied by 
£>Shen(z) Og, o.78/(£'Mill(z) Og. o.s), where Dshcn and Dmiii are the growth 
factor calculated for the cosmology used bv iShen et alj t2009,) and the Mil- 
lennium cosmology, respectively. 
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Figure 5. Excess bias between merger remnants and the parent galaxy pop- 
ulation (solid lines), for different stellar mass bins (the mass range is indi- 
cated in the upper-left corner of each panel in units of 10"'/i"'Mr:i). The 
dashed lines show the excess bias after matching the distribution of host 
subhalo masses. No excess bias is present if F = 1 (thin dotted line). 
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Figure 6. Median host DM subhalo mass corresponding to different stellar 
masses at different redshifts. The red bow-ties correspond to haloes hosting 
recently merged galaxies whereas blue triangles refer to the corresponding 
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distribution. 
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Figure 7. Bias of FOF halo merger remnants as a function of redshift for 
different mass ranges (as indicated on the plots, in units of IO '^/z^'Mq). 
The sy mbols indicate the bias of bright quasars c alculated by ISh en et alj 
120091) , inferred from the clustering observations o f IShen et al.l i2007ll 



tions iWvithe & Loebll2009h . The results of the previous sections 
for massive haloes and galaxies represent a challenge to this attrac- 
tive explanation, at least if the excess bias is to be of merger origin. 
Figure|7]suggests that quasars at z > 2 live in haloes ~ lO'^ /i^'Mq 
which is broadly consistent with the average host halo m ass esti- 
mated for lower redshift quasars (e.g.. lCroom et al.ll2005h . 

If bright quasars have no significant excess bias due to their 
merger-driven nature, as our results suggest, then either there is an- 
other unknown source of excess bias, or, more simply, their cluster- 
ing must trace the clustering of their host DM haloes and the dis- 
crepancy mentioned above must be explained in some other way. 
At this point it is therefore important to carry out a simple consis- 
tency check to see if there are enough massive haloes to host the 
luminous quasars actually observed in SDSS at z > 3. 

In Figure[8] we compare the number density of observed high- 
redshift quasars from Shen et al. (2007) with the number density 
of major halo mergers in the Millennium Simulation. Note that 
the information extracted from the simulation is a rate of merg- 
ers, i.e., the number of merger events within the time interval be- 
tween two snapshots (see Section |2j. Therefore, when comparing 
with quasar number densities we are forced to assume a quasar op- 
tical visibility time tq, that several independent studies have con- 
strained to be relativel y shor t and of the order of tq 10^ ~ 10^ 
yr (e.g.,j shankar et alIl20C)4 iMarconi et al.l |2004| ; iMartinil |2004 
lYu & Lull2004 IBirdet al.ll2008l , and references therein). In Fig- 



ure[8]we choose to multiply the rates by a quasar visibility time of 
lO^yr, which is, at those redshifts, the approximate time between 
two snapshots of the Millennium Simulation. The resulting cumula- 
tive number_densities are plotted in Figure[8]and are compared with 
the IShen et al] (l2007h quasar nu mber densities. The latter, taken 
from Table 5 in lShen et al. IIIoo^) and converted to our cosmology, 
are shown with a grey band in Figure [S] which takes into account 
a factor of two uncertainty due to possible sources not visible in 
optical surveys due to obscuration. 

From this plot we conclude find that if tq > 10** yr, there are 
potentially enough mergers of massive haloes in the Millennium 
Simulation to match the number density and the large-scale clus- 
tering properties of z > 3 quasars. A merger model would require 
a fraction 20-25% of haloes with mass > 8 X 10'^ /j-iM0 to be 
active at 3 < z < 4, in nice agreement with the analytical models 
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Figure 8. Number density of observed bright quasars from IShen et al] 
(20071) (gray line), compared with the number of major mergers in the Mil- 
lennium Simulation, obtained by multiplying the merger rate by a quasar 
lifetime tq = lO^yr. The solid blue lines refer to the cumulative number 
density of halo mergers, whereas the red dotted lines show the cumulative 
number of galaxy mergers. The minimum mass corresponding to each line 
is shown in the plot in units of [lO'^/i^'M;,] for the galaxies and in units 
of [10 *~/i~'M .,] for the haloes. The number densities quoted bv lShen et alj 
have been multiplied by a factor of two to account for objects miss- 
ing from optical surveys due to obscuration. 



of IShankar etd] j2008h . who find a duty cycle of 0.2 — 0.4 within 
the same redshift range. We stress here that our mapping between 
haloes and quasars neglects any scatter between halo mass and 
quasar luminosity, which could spoil the agreement as discussed by 
IWhite et"ail ( l2008h . Significant scatter in the Lqso — A^halo relation 
would decrease the bias of quasars, since many would be hosted by 
lower mass (and hence less clustered) haloes. 

The red-colored, dotted lines in Figure [8] mark instead the cu- 
mulative number densities of galaxy major mergers above differ- 
ent final masses, as labeled. The galaxy model predicts that, on 
average, the most massive galaxies that recently merged {Mq > 
8 X lO'" H^^Mq) reside in the most massive haloes of mass > 
8 X lO'^ /!^'M0, with a stellar-to-halo mass ratio consistent with 
the one empirically inferred from the cumulative number matching 
between the stellar and halo mass functions (e.g., fele & Ostriker 
2004; Sh ankar et alJl2006l : IConrov & Wechsleill20oi i^ lMoster etHI 
2009.) . However, we find that the number of major mergers for such 



massive galaxies is below the number of major mergers of the cor- 
responding hosts, as clearly seen in Figure[8]when comparing dot- 
ted to solid lines. 

This apparent discrepancy is explained, first of all, by the fact 
that when two haloes merge, their host galaxies will merge at some 
later time only if the new satellite halo loses enough mass to fall be- 
low the resolution limit of the simulation. Moreover, in the current 
treatment of galaxy mergers, when a galaxy becomes a satellite, it 
loses its hot gas component; cooling is then inhibited and the stel- 
lar component grows only moderately from the cold gas previously 
available. Therefore, although a given FOP halo merger may be 
counted as a major merger, by the time the corresponding galaxy 
merger occurs it may fall below our chosen threshold for a major 
merger. In any case, despite the fact that the number of major merg- 
ers is lower for galaxies than for haloes, the number of mergers of 
galaxies more massive than (Mp > 4 x lO'" A^'Mp) is still large 
enough to explain the observed number densities of bright quasars. 
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Figure 9. Upper panel: Bias of simulated bright quasars (B-band magnitude 
< —24 mag), compared with the bias of randomly selected subhaloes with 
the same mass distribution as the ones hosting the quasars. Lower panel: 
Excess clustering between the two populations. 



Taken at face value, the galaxy model would then predict that 
the SDSS luminous quasars detected at z > 2 should be hosted by 
galaxies as massive as > 4 x lO'** /j^'Mq. Given that virial rela- 
tions point to BHs more massive than > 3 x IO^Mq, this would 
suggest an increase, by a factor of > 3, of the BH-to-ste llar mass 
ratio with respect to local values iH aring & Rixl l2004h . In ad- 
dition, we find that the clustering of galaxies with stellar mass 
Mq ^ 4 X lO''* /i^'M© is too weak to match the observed quasar 
clustering. 

To better address the connection with the semianalytical 
galaxy models, we compare our results with the outputs of the de- 
tai led model for the co evolu tion of quasars and g alaxies presented 
in iMaruUi et al] ( |2008|) and iBonoU et al.l | |2009|) . Figure |9] shows 
the bias of luminou s optical quasars derived using the model of 
lMarumetaLl ( l2008l) . built on the a ssumption that quas ar activity is 
triggered during galaxy mergers. In lBonolietai] l l2009h we showed 
that such a model predicts well the clustering properties of observed 
optical quasars at a variety of redshifts and luminosities, indepen- 
dent of the specific light curve characterizing the active phase of a 
BH. In the upper panel of Figure|9l the bias of bright quasars in the 
model is compared with the bias of randomly selected dark-matter 
subhaloes with the same mass distribution as the ones hosting the 
quasars. The ratio between the two-point correlation functions of 
the two samples is shown in the lower panel: the excess bias is at 
most ~ 5%, except at z = 5, where the small number of simulated 
quasars results in a statistically unreliable result. 

It is clear that if bright quasars were hosted by DM sub- 
haloes less massive than inferred from the clustering analysis, the 
BH-to-stel lar mass ratio woul d be e ven higher (see also the dis- 
cussion in IShankar et all (e.g., |2008[) ). To address, in an indepen- 
dent way, the evolution of the average expected relation between 
the BH and its host, we compute the expected baryonic mass 
locke d in BHs follow i ng the method outlined by previous authors 
(e.g.. | Ferraresel I2OO2I; ICirasuolo et al.l [20051; IShankar et aklbood; 
IShankar & Mathuij|2007C IShankar et al.l2008n : first, we map haloes 
to their appropriate virial veloci ties V^\y at a given re dshift z apply- 
ing the virial theorem (e.g., Bar kana & Loebll200lh . We then link 
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by, e.g., iFerraresd ^20021) , and finally we compute the associate d 
BH mass via the local Mbh — d relation (e.g., Tun do et"al]|2007h . 
If we assume that these BHs are accreting at an Eddington ratio 
X > 0.5, we find that, at z = 4, all haloes above ~ 5 x IO'^/i^'Mq 
can indeed host a BH luminous enou gh to be recorded in the high- 
z quasar sample of Ishen et al.l ilOO^) . This simple exercise proves 
that if quasars are associated to normally biased haloes, the ratio 
between BH mass and halo mass could be similar to that observed 
locally. 



5 SUMMARY AND CONCLUSIONS 

In the present work we exploited the large halo and galaxy sam- 
ples extracted from the Millennium Simulation to test the idea that 
"merger bias", a tendency of recently merged systems to be more 
strongly clustered on large scales than typical systems of similar 
mass, could help reconcile the apparent discrepancy between the 
observed abundance and clustering of high redshift quasars with 
those predicted for massive dark haloes. Previous studies have, in 
fact, shown that the quasar number density and clustering can be 
simultaneously explained theoretically either by models character- 
ized by high duty cycles and negligible scatter in the Lqso — ^halo 
relation (White et al. 2008; Shankar et al. 2008), or by models with 
non-zero scatter and an excess bias for the haloes hosting quasars 
dWvithe & Lo"ebll2009l) . 

We quantify the importance of merger bias at different red- 
shifts and for different halo mass ranges. Defining as major merg- 
ers those events in which two friend-of-friend haloes of compara- 
ble mass merge into a single system between two simulation out- 
puts, we found that recently merged haloes with masses in the range 
5 X lO'^ to 1.6 X lO'^ /i^'Mp; show no significant excess clustering 
when compared to other haloes of similar mass. 

To connect with physically motivated models of galaxy forma- 
tion, we also looked for a possible merger bias among samples of 
galaxies selected from semi-an alytical models built on the Millen- 
nium Simulation JPe Lucia & Blaizot 2007). We considered galax- 
ies with stellar mass in the range 5x10^ — 1.6x 1O"/;^'M0 and 
found that merger remnants are typically 10 — 30% more clustered 
than other galaxies of the same mass. The merger remnants are 
hosted by systematically more massive subhaloes than other galax- 
ies, explaining a substantial part of this signal. However, even after 
correcting for this, we still observe excess bias at the level of ~ 5% 
for the most massive galaxy merger remnants, and of ~ 20% for 
our low-mass objects, which are insufficiently clustered to match 
the high bias of observed quasars. If we further restrict the analysis 
to central galaxies (i.e. galaxies at the center of a friend-of-friend 
group), for which a clear definition of halo mass is available, the 
excess clustering is once more diminished , approaching the null 
result obtained for haloes alone. 

The clear result obtained for haloes and massive galaxies in- 
dicates that merger bias is unlikely to be a viable solution to the 
apparent puzzle of the high clustering of high redshift quasars. On 
the other hand, we have also found that recently merged massive 
haloes withMhaio lO'^ /i^'Mq could be both numerous enough 
and clustered enough to match the observed quasar number den- 
sity and large-scale bias, if we assume a quasar visibility time 
?fy '--^ 1 X 10^ yr and if we assume negligible scatter in the relation 
between halo mass and quasar luminosity. 

In conclusion, if major mergers are responsible for triggering 
quasar activity, the lack of significant merger bias requires models 



to be characterized by high duty-cycles and negligible scatter in the 
relation between quasar luminosity and halo mass. 
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